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ABSTRACT 

We present a simple estimate of the mass "deficits" in cored spheroids, as a function 
of galaxy mass and radius within the galaxy. Previous attempts to measure such deficits de- 
pended on fitting some functional form to the profile at large radii extrapolating inwards; this 
is sensitive to the assumed functional form and does not allow for variation in nuclear pro- 
file shapes. For example, we show that literally interpreting the residuals from a single/cored 
Sersic function fit as implied "deficit" can be misleading. Instead, we take advantage of larger 
data sets to directly construct stellar mass profiles of observed systems and measure the stellar 
mass enclosed in a series of physical radii (M^,{< R)), for samples of cusp and core spheroids 
at the same stellar mass. We show that there is a significant (model-independent) bimodality 
in this distribution of central structure for this sample at small radii. We non-parametrically 
measure the median offset between core and cusp populations (the "deficit" AM^,{< R)). We 
can then construct the scoured mass profile as a function of radius, without reference to any 
assumed functional form. The mass deficit rises in power- law fashion (AM4< /?) cx/?''3-' **) 
from a significant but small mass at R< lOpc, to asymptote to a maximum ^ 0.5 — 2Mbh 
at ~ lOOpc, where Mbh is the mass of the central, supermassive black hole hosted by the 
spheroid. At larger radii there is no statistically significant separation between populations; 
the upper limit to the cumulative scoured mass at ^kpc is ^ 2 — 4Mbh- This does not depend 
strongly on stellar mass. The dispersion in M*(< R) appears larger in the core population, 
possibly reflecting the fact that core scouring increases the scatter in central profile shapes. 
We measure this broadening effect as a function of radius. The relatively low mass deficits 
inferred, and characteristic radii, are in good agreement with models of "scouring" from BH 
binary systems. 

Key words: galaxies: formation — galaxies: evolution — galaxies: active — quasars: general 
— cosmology: theory 



1 INTRODUCTION 

Massive ellipticals appear to exhibit central "cores" - regions of 
constant or weakly divergent surface brightness - as opposed to 
the power-law "cusps" - more ra pidly rising a t small radii - com- 
mon in les s massive elliptic als tond Il978l: lYoung et al.l Il978l: 
Lauei"l9 85l:lKormend'vlll98l.ll 987':'Lauer et alJll992l:ICrane et'al 



1993; La uer et al.l Il995l ; IFerrar ese et al. 1994 |2006| - ICote et al 



2006, 20 03). This bimodality in central s urface brightness slopes 
(Gebhardtet al.lll996l : iLauer et al.ll2007ij) participates in a larger, 
longer-recognized division of the elliptical population: the typi- 
cal giant, core elliptical rotates slowly, and has "boxy" isopho- 
tal shapes supported by anisotropic velocity dispersions; less 
massive cusp ellipticals (and SO galaxies and classical bulges) 
rotate more rapidly, have more isotropic velocity dispersions. 
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and exhibit "disky" isophotal shapes jKormendv & lUinaworthl 
'1982^; 'Davies et al.' 1983; Davis et al.' 'l985; Jedrzejewski et alj 


1987; Bender.. 1988.: 


Bender etal... 1989, 


1 992: .PeletieretaL.1 99(1: 


Kormendv & Bendei 1996: 


Faber et al.l 


19971: Simien & Prusniei 


il997a,b,c: Emsellem et al. 


2004, 2007 


: iMcDermid et al. 20061: 



e^ 



These differences thus led natur ally to the idea, devel oped i n 



Faber et al. ( 1997); Kormendv' (1999^; lOuillen et al 



Rest et al.l (12001.): ,Ravindranath et al., (.200lh: iLaine et al 



LaueretalJ j2005l 120073): IFerrarese et al.l ( |2006|) : ICote et al.l 
1 20071) : iKormendv et alj i^mi that disky, rapidly rotating cusp 



ellipticals are direct products of gas-rich ("wet") mergers, whereas 
boxy, slowly rotating core ellipticals have been shaped by sub- 
sequent dissipationless ("dry") re-mergers of two or more (ini- 
tially cuspy gas-rich merger remnant) ellipticals. However, nu- 
merical experiments and simple phase-space arguments demon- 
strate that central light profile shapes are, to lowest order, pre- 
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served in dissipadonless mergers without binary b lack holes (e.g. 
iBovlan-Kolchin et al.ll20o3 : iHopkins et Zll2009d) . It is therefore 
generally believed that the connection between the merger history 
of galaxies and their nuclear profile slope ari ses because of "scour- 
ing" by a binary black hole (fo r a review, see lGualandris & MerrittI 
l2007h . iBegelman et alj h980t) first pointed out that binary black 
holes coalescing in a dissipationless galaxy merger could stall (i.e. 
are no longer efficiently transported to the center via dynamical 
friction) at radii ~ pc, larger than the radii at which gravitational 
radiation can efficiently dissipate energy and merge the binary - 
the so-called "last parsec problem." They noted that significant gas 
content can provide a continuous source of drag and friction and 
solve this problem in gas-rich mergers, but that in "dry" mergers, 
the binary will remain stalled for some time, and will harden by 
scattering stars in the nucleus in three-body interactions. This will 
continue, flattening the nuclear slope, until a sufficient mass in stars 
(~ Mbh, by simple scaling arguments) is ejected to merge the bi- 
nary. 

It is therefore of particular interest to estimate the stellar mass 
which must be scattered to explain the slopes of cores, as a test of 
scouring models, a constraint on BH-BH merger rates and gravita- 
tional inspiral, and as a probe of the galaxy mer ger history. To date, 
there have been a few such calculations (see Milosavlievic et al.l 
2OO2I; iFerrarese et al.l[20o3: iLauer et al.ll20()7allbl: ICote et al JIiOOTI : 



Kormendv et a l.^'20Q9; Ko rmendv & Bendeill200 y): however, these 
have in some cases reached very different conclusions regarding 
the absolute mass "scoured," whether or not this scales with galaxy 
mass, and even whether or not there is any evidence for a bimodal 
population at all (as opposed to a smoothly scaling distribution of 
central profile gradients). 

Much of this ambiguity owes to the lack of an a priori physi- 
cal model for the "pre- scouring" profile shape. All of the estimates 
above implicitly assume some specific functional form for the light 
"scoured," where the parameters of this functional form are deter- 
mined from fitting the observed light profile at larger radii and ex- 
trapolating inwards. For exa mple, several of the above works fit ei- 
ther a "cored" Sersic model jOraham et al.l2003l) or a Sersic model 
with some deficit inside a "cutoff" radius to the large-scale galaxy 
profile, and then compare the actual light profile with the inwards 
extrapolation of the single Sersic portion of the fit (implicitly as- 
suming that the "pre-scouring" light profile followed the continua- 
tion of an identical Sersic profile to the outer, post-merger profile). 

However, there is no physical motivation for the Sersic model. 
It is well-known that there is considerable diversity in the cen- 
tral profile shapes within the cusp population - any individual 
system could deviate widely from this assumption - and more- 
over the assumption that profiles represent perfect Sersic con- 
tinuations is, in fact, known to break d own for the majority 
of non-core ellipticals (iFaber et alj Il997l: Iffibbard & YurJji99S : 



Rothberg & Joseplj2004lTruiillo et alJ2002l;lFerrarese et alj200e 
Cote et al ."2007: Kormendv et al."200g:'Hopki ns et al.ll2009bh . Foi 



For 

example, Hopkins et al. (2008b, 2009c) show that a physically mo- 
tivated decomposition into a two-component Sersic model more ac- 
curately approximates the central structure at ~kpc scales in ellip- 
ticals. Adopting a slightly different assumed functional form for the 
profile fit and extrapolation leads to large systematic differences in 
the "missing light" properties. In addition, it is well-established that 
core and cusp ellipticals have different profile shapes at radii out- 
side the "core" (see references above); even if all systems obeyed 
the assumed specific functional form, the use of the profile at large 
radii - radii that have no "information" regarding the scouring pro- 
cess - to anchor the fit could lead to significant differences between 



the true and inferred pre-sc ouring profiles (see lMilosavlievic et al.l 
l2002l : lHoDkins et al.l2009"dh . 

Moreover, the scouring process is not expected to be iden- 
tical from one system to the next. The efficiency and ef- 
fects of scouring, in detail, are functions of the merger his- 
tory (e.g. number of mergers and distribution o f merger mass 
ratios and orbits; see lOuinlan & Hemquisll 1 1991 iMerritll |2006| ; 
[Sesana et al. 2007h . as well as properties such as the triaxial- 
ity a nd stellar orbital anisotropy distribution as a function of ra- 
dius faollev -Bockelmann & Si gurdssonll2006l : lBerczik et alj2006l ; 
IPerets & A lexander 2008); numerical simulations suggest that it 
may not be described by a simple functional f orm relative to the 
profile at larger radii (see references above and IPerets et alj|2007l ; 
IzieJbOOTl) . 

In this paper, we therefore present a non-parametric estimator 
of the mass affected by BH scouring, as a function of radius. We 
construct stellar mass density profiles for a large sample of both 
cusp and core ellipticals, and use these to determine the stellar mass 
enclosed within different physical annuli. Considering this distribu- 
tion in ellipticals of similar mass, we show that - even without any 
reference to some "cusp" or "core" designation or profile fitting - 
there is a robust bimodality in the distribution of enclosed mass at 
small radii. We consider non-parametric measures of the median 
offset between these bimodal peaks; equivalently, the median dif- 
ference in enclosed mass between cusp and core ellipticals, as a 
function of radius. This allows us to construct - albeit still with 
some assumptions - the "scoured" mass profile, which grows with 
radius until asymptoting to a maximum scoured mass ~ 1 — 2Mbh 
near ~ lOOpc. We consider this as a function of galaxy mass, and 
with different methodologies, and show that it is robust. Unlike the 
approaches above, the methodologies here do not refer to specific 
profile fits, and account for the diversity of profile shapes within 
each population. Moreover this methodology allows us to actually 
measure the effects of scouring as a function of radius (rather than 
implicitly assuming them in the fitting); we show there is an ap- 
proximate power-law behavior at small radii with an asymptotic 
effect on the mass profile at larger radii. 

We adopt a = 0.3, = 0.7, Ho = TOkms"' Mpc"' cos- 
mology, but note that this has little effect on our conclusions. 



2 THE CONCERN: INFERENCE BY FITTED PROFILE 
SHAPES 

To begin, Figure [T] illustrates one reason why caution is needed in 
the approach of fitting galaxy profiles. Begin with a simple toy- 
model galaxy light profile, a Sersic profile with index = 3. This 
is typical of "cuspy," ~ L, ellipticals. Now imagine that the sys- 
tem undergoes some series of violent events, such that the profile 
at large radii is altered and an extended envelope is built up, but 
the profile at small radii is relatively unperturbed. This is a plau- 
sible scenario in the absence of core "scouring": in N-body and 
hydrodynamic simulations of "dry" mergers, even in major (1:1) 
mergers, the result to lowest order is that the profile shape in the 
central regions where the system is very dense tends to be con- 
served (stars there conserve their specific binding energies; see, 
e.g. Barnes 1988; Hernquist 1992, 1993), while some scattering or 
"splashing" of stars to larger radii makes the tail or envelope of the 
distribution at lar ge radii more extended (see lBovlan-Kolchin et al.l 
l2006l ; lHoDkinset al. 2009d). 

The effect is expected even more so if the envelope is built up 
not by a single very violent event, but by a series of minor mergers 
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Figure 1. Illustration of caveats in estimating the "missing" mass in core 
ellipticals via fitting some assumed functional form to the profile. Top: 
An idealized profile of a typical cuspy progenitor elliptical: here a perfect 
Hj = 3 Sersic profile, with the given effective radius. We plot surface stellar 
mass density versus radius. Middle: One possible way in which subsequent 
re-mergers and dry mergers will modify the profile, by building up an outer 
envelope of low-density material (here modeled with an n,, = 8 profile out- 
side of Rg). There is no change, however, in the central profile shape - in 
particular, no core "scouring" here. Bottom: The results of fitting these to 
assumed functional forms. First, a single Sersic law, fit down from large 
radii until the variance at a given radius becomes larger than the equivalent 
of = 0.1 mag. Comparing the inward extrapolation of this to the true 
profile gives an apparent "deficit" (mass fraction /deficit) inside of this ra- 
dius (Rcorc). Second, a core-Sersic profile fit (fitted to the entire profile). We 
plot the extrapolation of the outer Sersic portion of the fit (the best-fit full 
profile with "core" is indistinguishable from the true profile); with the given 
outer Sersic index, core break radius, and apparent core mass deficit. 



and/or mergers where the secondary is much less dense than the 
primary - in which case it is shredded at large radii and contributes 
directly to the envelo pe, but has a completely negligible effect on 
the central densities ([Gallagher & OstrikeJll972l : iHernguist et al.l 
ll993l : lNaabetal.l20oC Recent evidence from comparison of high- 
redshift massive spheroids and their descendants today suggests, 
in fact, that the buildup of such an envelope while central proper- 
ties are (relatively) weakly affected is probab ly the dominant evolu- 
tionary track for the most massive spheroid s jHopkins et alj2009al : 
iBezanson et ai]|2009l : iHopkins et al]|201(]h . We mock up such an 
expansion of the outer envelope by simply modifying the profile 
such that inside of the original R^, it follows exactly the same func- 
tional form as before; but at larger radii, it obeys a Sersic law with 
lis — 8 (these choices ensure the profile is smooth through the tran- 
sition). 

Now, fit the relic profile to a new Se rsic function. We consid er 
first a single Sersic function fit, following lKormendv et alj ( l2009h - 
we fit the profile at large radii first, and include more data points to 
smaller and smaller radii (re-fitting each time) until we reach some 
radius where the rms deviations about the inward extrapolation of 
the single Sersic law grow larger than 0.04 dex in log E, (roughly 
equivalent to A/i = 0.1 mag in surface brightness). Similarly, we 
can consider the core-Sersic profile 



E = : 



R% 



(1) 



which follows a Sersic function at large radius, with a break to 
some flatter power-law behavior at small radii (the fitted Rcok', 
Graham et al. 2003), fitted to all radii simultaneously. In either 
case, the resulting fit defines some outer Sersic profile, and some 
inner "core" region, inside of which the true profile falls below the 
inward extrapolation of the Sersic component. 

We show these extrapolations for both cases. They give simi- 
lar best-fit lis ~ 7 (obviously similar to the input ;is — 8 profile as 
r — >■ oo). The inwards extrapolation of each of these fits rises to 
much larger surface densities than the true profile. Thus there is an 
apparent large "core" or mass deficit inside small radius. Here, the 
apparent deficit appears at ~ 400 pc, and has a mass or light fraction 
(defined by simply integrating the light under each curve versus the 
extrapolated Sersic component) of ~ 2 — 5%. 

Of course, the inner regions of our mock profile are, by con- 
struction, identical to the progenitor - there has been no scour- 
ing of any kind. Assuming that the inwards extrapolation of the 
fit from large radii represents the progenitor is therefore mis- 
leading. Although the example here is purely illustrative, such 
an effect appears in numerical si mulations of dry galax y-galaxy 
mergers, for the reasons above. In lHopkins et al. I ( 12009 j) . the au- 
thors present a large suite of simulations of dry "re-mergers" of 
simulated elliptical galaxies, themselves the products of gas-rich 
mergers chosen specifically because of their very close agreement 
with the o bserved profile shapes of real observed cusp ellipticals 
([Hopkins et al. 2009b). After a 1:1 dry re-merger of such systems 
(with no scouring mechanism included), the authors note that many 
of the profiles exhibit exactly the sort of changes described above 
(where the inner profile is approximately conserved, but outer pro- 
file "puffed up"), and can be formally well-fit by a core-Sersic pro- 
file. This is despite the fact that there is no core scouring included 
or resolved in those simulations. Fitting each simulation with such 
a profile, we show the distribution of mass "deficits" (defined as 
above for the core-Sersic fit) and core break radii in Figure [2] The 
characteristic deficits of ~ 0.1 — 1% of Mgai, and break radii of 
~ 0.01 — 0.05 Re, are characteristic of the above issue. 



© 0000 RAS, MNRAS 000, 000-000 



4 Hopkins et al. 

1.01 



Simulations (No Scouring) 



0.8 - 



£ 0.6 - 



it 0.4 



0.2 - 



0.0 [ 



-1.5 



-1.0 



-0.5 0.0 0.5 
% Missing Light (fjefici,) 



1.0 



1.5 



1.0 r 



0.8 



^ 0.6 - 
o 

ra 

it 0.4- 



0.2 



0.001 



0.010 
Core Radius (R„ 



0.100 



e/Re) 



Figure 2. Distribution of "mass deficit" prope rties in 30 hydrodyn amic 
simulations of major dry merger remnants from Hop kins et alj (2009g). The 
galaxies are initially in good agreement with profiles of observed cusp el- 
lipticals, and are merged with no scouring mechanism. Post-merger, they 
are fitted to a core-Sersic profile, with the implied mass deficits and core 
break radii determined as in Figure [T] Top: Distribution of implied mass 
deficit percentages. Bottom: Distribution of imphed core break radii (Score 
in Equation[T] in units of Seff)- Owing to the effects in Figure[T] and simple 
variation/scatter in the true functional form of galaxy profile shapes, appar- 
ent non-trivial mass deficits and large break radii can arise. 



The problem arises because of the nature of the Sersic pro- 
file, but could generally arise from any extrapolation from larger 
radii. Fitting such a profile implicitly couples the profile shape at 
large radii to that at small radii. For a Sersic profile, a more ex- 
tended envelope at large radii corresponds to a larger index Wj. 
However, a larger n, also corresponds to a more steeply rising pro- 
file at small radii, as is clearly evident in Figure [T] The effects 
on the inferred nuclear "scoured" mass can be large: an iij = 2 
profile has just 0.06% (3%) of its mass within 1% (10%) of Re, 
an ns = 8 profile has 2% (14%). Thus, any change in the profile 
shape at large radii will necessarily change the inferred "progeni- 
tor" profile at small radii. Moreover, observed core ellipticals tend 
to have significantly higher iis in their outer reg ions (more extended 
envelopes) than lower-mass cusp elliptical s dTruiillo et al.ll200^ : 
iFerrarese et alj|2006l : iKormendv et alj|2009l) . This means that, as- 
suming such an inward extrapolation, the assumed "progenitor" 
small-scale profile in the cored systems is in fact more steep than 
that observed in any actual cusp elliptical. For exa mple, the stel- 
lar mas s profiles of simulated merger remnants in iHopkins et al.l 
( l2009bl. [2008a) are not generically reproduced by a single Sersic 
function, but rather exhibit multi-component structure indicative of 
the mixed roles of dissipation and violent relaxation. 



3 METHODOLOGY: AN ESTIMATOR OF THE CORE 
MASS "DEFICIT" 

Given these concerns, we desire here to develop a more simplified 
but hopefully more robust check of mass deficits and core radii. 

We begin by considering the surface stellar mass density pro- 
files as a function of projected radius R for ellipticals and 
SO galaxies with HST imaging of their nuclei (and ground-based 
profile data at large radii, allowing accurate surface brightness pro- 
file measurements from < lOpc to ~ 50kpc). We compile ob- 
served surface bri ghtness profiles from lKormendv et al.l ( |2009[) and 
iLauer et^ ( l2007ai) : this consists of a total of ~ 180 unique lo- 
cal ellipticalsQ The HST data allows robust classification of the 
cusp/core status in each system, discussed therein. We determine 
stellar masses for each system in , Hopkins et al., (2009b), based on 
the integrated observed optical luminosity in several bands, us- 
ing the m ass-to-light ratio cal ibrations as a f unction of color in 
iBell et al.l |2003) corrected to a lChabrieJ j2003h IMfQ The isopho- 
tally averaged major axis profiles are measured in rest-frame opti- 
cal; we convert to a stellar mass profile based on the measured total 
stellar masses and the assumption of a radius-independent stellar 
mass-to-light ratio. Most of the objects also have resolved color 
gradients - we find that converting to a stellar mass profile using the 
local color and a color-dependent M* /L makes no difference, as the 
color gradients are weak. Likewise, conversion to stellar mass pro- 
files using stellar population gradients and comparison of profiles 
from different instruments and wavebands in these samples are dis- 
cussed extensively in Hopkins et al. ( 2009b, 2008a); the differences 
are much smaller than the variation between individual profiles, and 
do not affect our conclusions. The .Kormendv et aL 1 ,200ft) sample 
is a volume-limited survey of the Virgo spheroid population; as 
suc h it incl udes few very massive galaxies (M* > 3 x lO" Mq). 
The lLaueT et al. ( 2007a) galaxies are chosen to be characteristic of 
massive ellipticals in the local Universe, including more massive 
systems up to a couple lO'^ Mq . At the masses of interest, both are 
representative of the distribution of spheroid sizes, velocity disper- 
sions, and profile shapes in the local SDSS galaxy sample (see e.g. 
iHyde & Bernardi 2008). 

Figure [3] illustrates the methodology. First, consider a narrow 
range in total galaxy stellar mass (here M, = 10""^"'' Mq), in 
which there are both cusp and core elliptical populations. Ideally, 
we would also consider a specific range in black hole mass, but di- 
rect BH mass measurements are available for only a small subset 
of these objects (too few for the statistics here); as s uch, we adopt 
the observed BH-bulge stellar mass relation from i Haring & Rixl 
(2004) to estimate the corresponding BH masses in these hosts 
(more on this below). We show a direct comparison of the sur- 
face stellar mass density profiles of these systems. We separately 
denote the cusp and core ellipticals|3 We plot results down to the 
minimum quoted radii from the sources above (the minimum radii 



' Note that although the c omposite (HST+ground-based) profiles were 
used in lLaueretaiTj2007iJ) to estimate effective radii, they were not ac- 
tually shown in the paper. 

^ Varying the specific bands used to determine stellar masses makes little 
difference, and changing the IMF will systematically alter the stellar masses 
of all objects considered, but will not affect our comparisons. 
^ For co nsistency, we adopt the c usp/ core designa t ions in the source sam- 
ples from ' Kormendv et al J )2009h and lLauer et li] j2Q07 j) . based on their 
logarithmic nuclear profile slopes at the smallest radii observed. However, 
other classifications are identical for all but a few objects; adopting those in 
iFerrarese et aU l ,2006i) or .Cote et al., (2007) makes no difference. 
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Figure 3. Top Left: Stellar mass surface density profiles of spheroids in a narrow stellar mass range Mt = 10"-" ^Mq. Cusp and core ellipticals are 
separately noted. Top Right: Histogram of the enclosed mass M« {< R) inside some radius, determined from the surface density profiles, at S = lOpc. There 
is a clear bimodality. Gaussian fits to the cusp and core populations separately are shown (histogram of each in dotted lines). Bottom Left: Same, at 1 kpc. 
Bottom Right: Difference AM* between the median enclosed mass M»(< R) at each R, comparing the cusp and core populations (black circles). Points 
represent detection of a statistically significant offset in the median M, (< R) between the populations, with the absolute mass noted. Error bars extending 
below the plotted limits should be considered upper limits; in this case the observations are consistent with no difference. Magenta squares fit the histogram to 
a double-Gaussian, i.e. assume no prior cusp/core designation (as such the uncertainties are larger); they are plotted only where a double Gaussian is favored 
at any significance over a single Gaussian. Compare the average BH mass in galaxies of this stellar mass range. There is a clear bimodality at small scales, 
where AM» (< R) grows with R until appearing to asymptote to a maximum ~ Mbh at R ~ 100 pc. 



at which their seeing-deconvolved profiles are considered reUable) 
- essentially the HST resolution limit. We stress that we will use 
these points in all of our analysis below, so some caution is needed 
in considering how robust the conclusions we quote below are at 
the smallest radii. In general, resolution limits are such that below 
~ lOpc, the sample size with any resolved points drops dramati- 
cally. In the more massive galaxies, the limits are larger; the most 
massive systems generally are limited to > 30 — 50 pc resolution. If 
we simply interpolate the nuclear profiles inwards based on their 
slope at the minimum resolved radii, our results remain similar 
down to arbitrarily small radius, but we are specifically attempt- 
ing to avoid such assumptions about the "true" profile shape, so in 
general one should consider the smallest-radii points in each mass 
bin to be fairly uncertain. 

Bearing those caveats in mind, at some fixed radius R, we can 
consider the distribution of stellar mass enclosed within that radius, 
i.e. 

nR 

M,(<R)= E{R')2TvR'dR' . (2) 
Jo 



Note that this is for the projected surface density E. We can also 
apply an inverse Abel transform to de-convolve the (circularized) 
projected density profile to recover the three-dimensional density, 

1 f"" dE dR 

P(r) = - - / , (3) 

' TV J, dR 

dSracewelill 1 965h and then integrate to determine M, (< r) as the 
mass within some three-dimensional spherical radius. Because M« 
is an integral quantity, this makes little difference; in what follows 
we prefer the use of M, (< R) as it is more numerically stable (the 
gain in accuracy in de-projecting to a three-dimensional profile 
is offset by the fact that the Abel transform required in that de- 
projection is very sensitive to small uncertainties in the measured 
profile; as such a careful error analysis shows no improvement in 
the constraints). Figure |3| shows this for the systems in the given 
total mass range, at ^ = lOpc and R = I kpc. 

We first consider the cusp and core systems separately; Fig- 
ure|3]plots the best-fit Gaussian to each distribution (i.e. Gaussian 
with the same median and 1 a dispersion), and the cumulative dis- 
tribution (i.e. histogram of all sources, without reference to their 
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cusp/core status). At small radii ^ lOpc, there is a clear difference 
between the distribution of cusp and core enclosed masses - signif- 
icant at ~ 5 (7. This is clear directly in the total histogram - there 
is a striking bimodality in the enclosed mass. Comparing a double 
versus single normal distribution, we find that a distribution with 
two distinct peaks is preferred at > 4a. Other tests yield similar 
results; for example, the "dip" test of 'Hartigan & HartigaiJ ( flisl) 
favors bimodality at ~ 3 — 4 cr without reference to an assumption 
of normality in the two-peaked distribution, and can be used to non- 
parametrically quantify (albeit with larger uncertainties) the peak 
separation, yielding similar results. This should not be surprising; 
the statistical significance of the bimodality in central galaxy den- 
sitie s (and their slopes) has been the subject of a number of works 
(see 



> (and their slopes) has been the sub|ect or a number ot works 
Ferrarese etal]|l994 lOebhardt et a"l]|l996l ; iFaber et alJ[T997l : 



iLauer et alj|2007al) . and is well established in at least this sample 
of chosen cusp-core galaxies; whether it is smoothed out in strict 
volume-limited samples (the degree of "dichotomy" or the unifor- 
mity of the trans ition) is however a question still being debated (see 
ICote et al.ll2007h . On the other hand, however, at 1 kpc, the distri- 
bution is clearly unimodal, and there is no significant difference 
between cusp and core profiles. 

Given this bimodality, we can quantify the difference, on av- 
erage, between the mass enclosed at small R in the populations. 
Specifically, we can take the absolute difference between the me- 
dian M, (< R) for cusp and core ellipticals, with the appropriate 
uncertainty reflecting both the uncertainty in each median and the 
weighted uncertainty in the difference. This is non-parametric, and 
allows us to directly compare cusp and core profiles at a given mass 
without reference to some model for what the profile "should be." 
Moreover, the large scatter in M, (< R), reflecting the scatter in 
profile shapes, is properly accounted for - we can robustly quan- 
tify the median mass difference in the populations. And it also has 
the advantage that we can robustly define an error bar on the mass 
"deficit." 

Figure [3] shows the result; on this scale, the two populations 
are offset by a median AM.(< R = lOpc) = l.lloj x 10** Mq. 
Compare this to the average BH mass expected and observed in 
spheroids of this total stellar mass range, (Mbh) ~ 2.5 x 10^ Mq 
(using Mbh ~ 0.00 14 M., from Haring & Rix 2004). 

We can be even more conservative, and make no refer- 
ence to the prior designation of systems as "cusp" or "core." In- 
stead, we simply fit the cumulative distribution of M, (< R) to 
a double-peaked (double Gaussian) distribution (equivalently, we 
non-parametrically measure the double-peak separation from other 
statistical tests as discussed above). In this case, the "separation" is 
only well-defined if the double-peaked fit provides a significantly 
better fit than a single-peaked fit; if a double-Gaussian fit is signifi- 
cantly better (with the attendant 6 degrees of freedom instead of 3), 
i.e. results in both Ax^ > 1 and improved /v, then we quantify 
the resulting difference in the median between the two fitted peaks 
and the formal statistical uncertainty associated. Unfortunately, be- 
cause the number of objects in any mass bin is small, we find that 
freeing all the parameters of the fit yields too many degrees of free- 
dom for any but marginally significant results (although the results, 
averaging over all mass bins, reach ~ 3 a significance and support 
our conclusions from the prior-constrained comparison of cusp and 
core populations above). 

We therefore show results fixing the relative fraction of the 
population in each Gaussian to the fraction of cusp or core ellip- 
ticals (we have also fixed the relative widths, rather than the frac- 
tional population, and find similar results). If the division between 
cusp and core populations did not reflect some physical bimodal- 



ity, this measure would give a null result (no signal for statistically 
significant difference) at all radii; if a large fraction of objects were 
somehow misclassified, or if an "intermediate" class exists (i.e. a 
smooth transition from cusp to core), this method will allow for 
that (even with our fixing of either the Gaussian widths or relative 
normalizations), and should still robustly quantify the difference 
between bimodal peaks in the population. We find similar results to 
the case above, but with larger uncertainties. 

Figure [3] continues to show these fit results as a function of 
radius, from ~ 3pc to ~ lOkpc. Again, we recall the cautions 
stated above regarding the effects of resolution and seeing limits 
at the smallest radii. At > lOOpc, however, the distribution of en- 
closed masses is clearly unimodal. Specifically, at ~ 1 kpc (where 
we show the full distribution) for example, fitting directly the cu- 
mulative distribution there is no preference for a doubly-peaked 
underlying distribution^ Even considering the cusp and core popu- 
lations separately, there is no statistically significant difference be- 
tween their M*(< R) distributions on this scale. Formally, we find 
AM. = -3 X 10' Mq (+1.7 X IO' Mq) (-2.0 X 10** Mq); in other 
words, the distributions favor no difference between the two pop- 
ulations, and we can set a 90% upper limit on the core "deficit" at 
this radius of ~ 4Mbh. 

In Figure|4l we generalize this to several mass intervals, from 
~ 10" - 10'^ Mq. In each of several mass intervals, we again 
show the mass profiles of the observed cusp and core systems at 
small radii. We also show the logarithmic slope of the mass pro- 
file, which highlights the flattening in the cusp systems (discussed 
further below). And we show the median AM, (< R) as a func- 
tion of radius R, determined using the same methodology as in 
Figure [3] To put the different mass intervals on the same foot- 
ing, we show the mass fraction M, (< R) /M «(< R — » oo). In these 
units, the average BH mass implied by the iHaring & Ri}J i2004h 
BH-host galaxy mass correlation is a constant ~ 0.0014. That is 
not, of course, the only determination of th is relation (compare 
Ma gorrian et al ] |l998l ; iMarconi & Hun3l2003h , nor is it the unique 
predictor of BH mass: ve locity dispersion jFerrarese & Merrit^ 
l200d ; iGebhardt et al.ll2000h. host gala x y binding energy or poten - 
tial well depth /Hopkins et al."2007lJH lAller & Richstonell2007h . 
and profile shape ( Graham etal. 200 1|) are possible alternatives. 
However, given the relatively limited mass range shown, these pre- 
dictors do not differ very dramatically over this range. Specifically, 
we also show the median Mbh /M« expected in each bin using the 
Mbh — cr relation together with the measured cr of each galaxy in 
the subsample. At the highest masses, the median expectation from 
Mbh — (7 is a factor of ~ 2 lower than that from the Mbh — Mbuige 
relation, but this is where our error bars become quite large (owing 
to the dearth of cusp ellipticals at these masses), so our comparison 
cannot discriminate at any significance whether there is a better cor- 
relation between the implied deficits and the BH-host mass implied 
BH masses or BH-u implied masses. It should however be bom in 
mind that the two estimators will diverge for the most massive sys- 
tems, near the upper limit of o ur sample (by M, = lO'^ Mq they 
may differ by a factor of several. lLauer et alj2007bl . but most of the 

At these radii, caution is needed adopting the fully freed bimodal fitting 
approach - if the distribution is slightly non-Gaussian (has some skewness), 
and one fits a double Gaussian, it is possible to obtain a statistically signif- 
icant "offset" that is clearly unrelated to the cusp/core populations. Some 
constraint like that adopted here is necessary, or more arbitrary functional 
forms (still unimodal) must be fit in combination, to prevent unphysical 
results. Given these uncertainties, properly marginalizing over the fitting 
parameters is particularly important to obtaining valid error estimates. 
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Figure 4. As Figure[3] for spheroids in different stellar mass intervals. Top: Mass profiles. Middle: Logarithmic slope of the mass profile, which highlights 
the degree of "flattening" in the central regions. Bottom: Estimated "deficit" as in Figure |3] as a function of radius. To compare systems with different total 
masses, we plot the fractional mass difference enclosed, i.e. AMt(< S)/M« (normalizing each profile by its stellar mass; because of this, the plotted results 
for the ~ 10"~"''Mq bin are not exactly identical to those i n Figure[3] though the y are very similar). Dotted lines show the expected (Mbh) for the median 
galaxy stellar mass in the bin and the BH-host mass relation of lHaring & Rixl )2004l) . or the median velocity dispersion a and BH-cr relation of lTremaine et alj 
(200Z) (blue and green, respectively). In each mass bin, similar trends are seen: a significant bimodality at small radii with AA/, (< R) growing to a maximum 
detection at ~ 0.5 — 2A/bh at S ~ lOOpc; in at least several mass bins, 90% upper limits of AM» < 2 — 4Mbh are present out to ~ 0.5 kpc. 



objects in our most massive bin are concentrated near the low-mass 
range of the interval). 

In each case in Figure|4l a similar trend is seen. At very small 
radii < lOpc, there is a significant bimodality in the population, but 
the amount of mass contained in this radius is small. As such, the 
absolute mass difference between the cusp/core populations grows 
to larger and larger radii (very roughly with a power law-index 
AM* oc R^'^^^'^). However, between ~ 10— 100 pc, this appears to 
asymptote to a maximum of AM, ~ 1 — 2 Mbh- Comparing each 
panel with the appropriate implied (Mbh) from either the BH-host 
mass or BH-tj relation, AM, appears similar in units of M, . At 
larger radii R > lOOpc, there is no statistically significant bimodal- 
ity remaining. Formally, the difference AM, — >■ at these radii; if 
this were strictly the case, it would imply that scoured mass is sim- 
ply "kicked" to slightly larger radii (not expelled from the galaxy 
or moved to very large radii). However, the error bars are still large 
- it is very possible that the mass is removed entirely, and so the 
mass difference AM, remains constant at ~a couple Mbh at all 
R (of course, this will become less significant relative to the total 
enclosed mass at larger R). But the difference does not appear to 
grow; the ~ 90% upper limit in at least some cases remains at or 
below ~ 4Mbh at radii ~ 0.5 — 2kpc. 

Of course, the comparison of profile shapes makes it immedi- 
ately clear that the uncertainties at both low and high masses are 
considerable, owing to the lack of core and cusp spheroids, respec- 
tively. Note that the total mass range in Figure|4]is relatively small, 
and the individual bins in mass are relatively wide. This is neces- 
sary for the statistics needed, and because cusp galaxies become 
progressively more rare at high masses as cores are more rare at 
low masses. But together with the limited sample sizes observed 
at sufficiently high resolution, it means that the difference between 
the median mass of the cored spheroids in our lowest-mass bin and 



that of the cusp spheroids in our highest-mass bin is only a fac- 
tor of ~ 7. And recall that the apparent cusps or cores at the mass 
"extremes" may also represent unusual growth histories rather than 
typical progenitors/descendants. The effect of resolution limits is 
also increasingly evident in the highest-mass interval, where the 
innermost point of any comparison is at ~ 20 — 30pc (with the best 
statistics at ~ 50 — lOOpc). Thus our results at the extremes should 
be taken with appropriate caution. 

Figure|4]also shows the logarithmic derivative of the mass pro- 
file, — rflogE,/dlog/?, as a function of radius. This highlights the 
range of behavior in the central regions of the observed systems. It 
also reaffirms a number of our results. At large radii, it is clear that 
it is not just the absolute values of E, that are similar in all systems 
in each bin, but also their slopes - this is where we do not see any 
statistically significant separation. At small radii, however, we see 
what is expected - the galaxies identified as "cusp" systems gener- 
ally maintain slopes in the range —1 < dlog E,/<ilogi? < —0.5, 
down to the resolution limits, whereas those identified as cores 
fall closer to flat slopes. Of course, there are also intermediate 
systems in between. IVlost interesting, the radii at which we see 
AM, appear to asymptote (or above which do not see a statisti- 
cally significant difference) appear to roughly correspond to the 
radii above which the core populations reach slopes steeper than 
d\ogT.-M / dlogR < —0.5. In the most massive core systems (al- 
though it is more ambiguous in the less massive core systems) this 
is also typically within a factor ~ 2 of the radius where the abso- 
lute value of the second derivative of the profile slope is maximized 
(if it is estimated with some smoothing kernel; we do not show the 
numerical second derivatives here because they are very noisy). 

This is interesting because such a threshold has been proposed 
and used as an estimator of where "scouring" action has taken place 
- if valid, it would allow an assessment of the scouring effects in in- 
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Figure 5. Difference in the logarithmic dispersion in M»(< R) in cusp 
and core populations ((7[log{M, (< R)}], for each population, subtracted 
in quadrature). Where this is > 0, it implies a significant broadening of the 
dispersion in M«(< R) by the plotted dispersion. Style is as Figure |4] we 
exclude the two low-mass bins as they show no statistically significant de- 
tection at any radius. We find no statistically significant cases where the 
dispersion in core systems is more narrow than that in cusp systems. The 
observations tentatively show evidence for non-trivial variation in scouring 
processes at small radii. 



larger for the core population, as compared to the cusp population. 
In contrast, there are no cases of statistical significance where the 
cusp population has larger dispersion. This might be anticipated, 
given the arguments in §[T] the scouring process is, in detail, sen- 
sitive to quantities such as the merger mass ratio, binary orbital pa- 
rameters, and triaxiality of the host, that will vary even at fixed stel- 
lar mass and fixed host mass profile shape. As such, it is plausible 
that scouring will broaden the distribution of mass profile shapes 
and correspondingly M,{<R) at small radii, relative to the distribu- 
tion present in the progenitor population. If both a [log {M, (< R)}] 
for core and cusp populations can be determined, this additional 
broadening or dispersion term should be given by their difference 
(subtracted in quadrature). Figure [5] shows this. Because this re- 
quires going to higher-order moments than the median Mt(< R), 
the statistics are correspondingly more limited - we only find a sta- 
tistically significant result at a couple of radii and masses. But in 
these cases, the trend is interesting. Reflecting what is seen in the 
absolute mass differences, the "additional" dispersion term appears 
to grow at small radii (although it should be stressed that the disper- 
sion in both populations grows at small radius; see lHopkins et al.l 
[2009 e). At ~ 1 — lOpc, scattering processes may lead to fractional 
dispersion in AM, (< R) at the ~ 0.5 dex level. In absolute terms, 
this amounts to dispersion, at these radii, of ~ Mbh in the absolute 
mass differences. At larger radii, the fractional dispersion falls (as 
it must, if the absolute mass scattered does not significantly exceed 
~ Mbh). 



dividual galaxies, without the l imitation of large s amples required 
here fsee lOebhardt et al ] | 19961 : iLauer et al.ll2007bh . Of course, we 
are expressly attempting to avoid making assumptions about the 
functional form of the pre-scouring profile, so a comprehensive 
comparison of these results to those here is outside the scope of this 
paper, but the comparison here may lend some motivation to those 
approaches, if the threshold slope values are taken from compari- 
son similar to that above. More rigorously, we can ask: if we were 
to adopt some parametric extrapolation of the profile shape on a 
per-galaxy basis, what would provide the best match to our results? 
Consider each mass interval above. For each core system, we could 
assume (with the dangers outlined above inherent in this process) 
that the pre-scouring profile followed a power-law at small radii 
with approximately constant slope at some dlog 'E^/dlogR — rjaiu 
we simply would estimate the "pre-scouring" profile by extrapolat- 
ing the observed profile with rjcnt at radii below that where the ob- 
served profile slope falls below r^crit. Averaged over the core popula- 
tion in each mass interval, this will give some average AM, (< R) 
implied ~ and we can ask what value of 77crit minimizes the dif- 
ference between this estimator and our non-parametric estimator. 
Performing this exercise for the mass intervals in Figure |4l we 
find that in all mass bins, the difference between our estimator and 
the distribution of deficits recovered by this method is minimized 
for threshold slopes in the range —0.5 < r^crit ^ —0.8 (specifically 
—0.75, —0.52, —0.81, —0.63, —0.60, from least to most massive, 
with respect to the cusp vs. core estimator; the error bars are suffi- 
ciently large in the 2-Gaussian estimator that there is not a strong 
discriminant in this range). Considering all bins simultaneously, the 
best-fit value is rjait ~ —0.62, for this sample. 

At each radius, we have not just the median M, (< R) for each 
of the cusp and core populations, but also the dispersion (roughly 
lognormal) in M, (< iJ), a [log {M, (< i?)}]. At several radii in the 
mass bins considered, this dispersion appears to be significantly 



4 DISCUSSION 

We illustrate important systematic uncertainties in the mechanism 
of inferring the "mass deficits" at the centers of cored spheroids 
via fitting of prior assumed functional forms to the mass profile 
(the "core" in the "core-Sersic" profile), and show that these can 
in principle lead to large artificial mass deficits and large apparent 
core radii. Specifically, we show that the assumption that the "un- 
scoured" progenitor profile matches the inwards extrapolation of a 
single Sersic fit to the galaxy profile at larger radii can be prob- 
lematic, and is often misleading when applied to simulated galaxy 
profiles. 

To avoid these uncertainties, we therefore define a simple, 
non-parametric measure to compare the stellar mass difference in 
the center of cored versus cusp spheroids. Essentially, this amounts 
to a comparison of the enclosed mass profiles M, (< R), quanti- 
fying the significance of bimodality in this distribution and, given 
that, the separation between bimodal peaks. This does not rely on 
any prior assumption regarding the functional form of the "true" 
pre-scouring profile (only the concept that core ellipticals descend 
from cusp ellipticals similar to those present today); nor does the 
inferred mass deficit depend on the profile behavior at larger radii. 
It also allows us to account for the non-trivial dispersion in central 
profile shapes within the cusp and core populations. 

We stress, however, that we are still making a very impor- 
tant assumption: namely, that the cusp ellipticals present today, in 
a given stellar mass bin, are representative of the true progenitors 
of the core ellipticals present today at similar stellar mass. This is 
not necessarily the case! The progenitor galaxies had to form earlier 
than their descendants (and the stellar populations of the core popu- 
lation appear older, on average, than those of the cusp population). 
So it is not hard to imagine, for example, that the core progeni- 
tors formed from most gas-rich events at high redshift, and/or that 
they were more concentrated and compact, as observed spheroids 
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at these redshifts (e.g. lWilliams et alj201(]| . and references therein). 
And some of the cusp population at high masses may represent 
cored galaxies that had their cusps "rebuilt" by new inflows from 
subsequent gas-rich mergers and/or cooling flows. Indeed, struc- 
tural evolution in the outer parts of these galaxies, and subsequent 
merg e rs, are expected in most theoretical models jHopkins et al.l 
l20ld iNaab et alj |2009^ It does appear, however, that the central 
portions from ~ 1 — 5 kpc of even the most compact z > 2 spheroids 
are similar in their stellar mass densities and profile shapes to local 
massive spheroids ("Hopkins et alj |2009ah . so our implicit assump- 
tion may be justified. Nevertheless, an ideal comparison would 
contrast present-day cored spheroids with higher-redshift candidate 
progenitors, chosen on the basis of a combination of mass and clus- 
tering or other constraints such that they are likely to evolve into 
the present-day core population. Unfortunately, there are no present 
observations that can resolve the cusp/core structure of these high- 
redshift systems. Therefore we are limited to the present day cusp 
population as our best-guess reference set, and should be careful to 
interpret our results in this context. 

With that caveat in mind, we show that there is a clear trend in 
the comparison of present-day nuclear profiles: at small radii, the 
distribution of enclosed mass in chosen cusp-core samples is bi- 
modal at high significance. The mean mass difference - the average 
"scoured" mass, in the scouring interpretation - between the two bi- 
modal peaks is low at small radii Mbh at ~ 3 — lOpc), reflect- 
ing the fact that there is simply little absolute mass in such a small 
enclosed volume. The apparent "scoured" mass then grows with 
radius (roughly as cx i?'^"' *) until ~ 10- lOOpc, at which point 
it asymptotes to a constant mass fraction ~ 1 — 2Mbh (~ 10~^^ — 
10^^ 'M,). The fact that we see this "turnover" or asymptotic be- 
havior suggests that we are actually resolving the total scoured 
mass, and radii affected. At larger radii 2> 100 pc, there is no statis- 
tically significant signal for bimodality, and the ~ 90% upper limits 
on the scoured mass remain consistent with the same asymptotic 
mass difference. Neither the scoured mass fraction nor characteris- 
tic radii appear to scale strongly with galaxy mass; however, our dy- 
namic range is limited to the mass range over which both cusp and 
core populations exist in sig nificant numbers (~ 10"'^~"^Mq). 
The characteristic scouring radius, for example, may well scale 
with the BH radius of influence Rbh ^ GMbho""^ cx M°'^, and 
in Figures I4I5I there is a tentative hint of this (it appears that the 
characteristic radius increases with mass), but the difference is not 
statistically significant over this dynamic range with the present 
statistics. 

These results are consistent with the expectation from A'-body 
experiments of the effects of a binary BH on the central stellar 
mass distribution. Such experiments typically find that the stellar 
mass scattered before the binary BH merges is ~ 0.5 — 1.5 Mbh 
(with dependence in that range on the central structural proper- 
ties of the galaxy and structure o f the inspiralling mass; see lMerrit^ 
l2006l ; |Perets & AlexandeJl2008l and references therein). The char- 
acteristic radii are also consistent; such simulations typically find 
that although stars from larger radii can contribute significantly 
to BH coalesce nce i f they are on long rad ial box orbits in a tri- 
axial pote ntial jOuin lan & He rnquislll 19971) or if there is a strong 
perturber jPerets et aU i200 7). a significant impact on the mass 
profi le is restri c ted to radii where the initial M, (< R) < Mbh; 
for a iHernquistI fT990) or r ' bulge, this is approximately below 
R/Rcf! « 0.5 \/ Mbh /Mast or ^ lOOpc for an ~ L, bulge. 

Comparing this to studies where the mass "deficit" AM, 
is determined from profile fitting, w e find a smalle r AM, (and 
smaller scoured radii) than some. iFerrarese et al. I i2006h and 



iKormendv & Bended ^20091) estimate AM, in a sub-set of the 
galaxies here, by comparing the actual light profile to the inward 
extrapolation of Sersic functions fitted to larger radii, they find 
AM, ~ 2 — 4 Mbh and ~ (10 — 20) Mbh, respectively. However, 
we show in §|2]that these differences can arise as a systematic ef- 
fect of the assumed profile shape. Fundamentally, at present, any 
such choice is somewhat arbitrary. And by assuming a specific 
functional form, the fits implicitly couple the assumed "progeni- 
tor" profile at small radii to the profile at large radii; for the Sersic 
fit, a galaxy with a more extended envelope (higher n,) will nec- 
essarily imply a more steep nuclear Sersic profile extrapolation. 
Thus, for the same true nuclear profile, the inferred mass deficit 
and core radius increase. And in fact, in both samples, the observed 
core ellipticals have more extended envelopes and higher than 
essentially any of the cusp elliptical population - thus the inward 
extrapolation of their Sersic profiles is not necessarily identical to 

the t ypical nucle a r shape in c usps. 

iLauer et alj ( l2007bl) and iMilosavlievic et alj ( l2002h consider 
instead the absolute mass enclosed within the radius where the pro- 
file steepens beyond a specified logarithmic slope (dln/^/dlnr = 
—0.5) or relative to the inwards extrapolation of a power-law (as 
opposed to Sersic) fit. They find AM, « 2.5Mbh and ~ 0.5 — 
IOMbh, respectively. Both cases, however, still require some im- 
plicit a priori assumption of a "correct" functional form for the 
progenitor light profile at small radii (or the characteristic shape of 
"scoured" regions). Moreover, in this approach, the mass measured 
is really the post-scouring mass inside this radius, not necessarily 
the mass removed (a non-trivial difference at the factor ~ 2 — 4 
level). 

In addition, the nature of the relatively flat profile at small radii 
means that whenever some functional form for the progenitor light 
profile is assumed, most of the "missing light" will come from the 
radii very near where the core breaks downwards relative to the 
"non-core" profile. As such, the inferred mass deficits are quite 
sensitive to the exact parameterization of the extrapolated profile 
near this radius (or, for a fixed slope cutoff, to the exact slope cut 
chosen). For the core-Sersic law, for example, there is a degeneracy 
between the fitted core break steepness parameter and the inner pro- 
file slope; small variations and common choices in fixing or free- 
ing one of those parameters systematically affect the inferred mass 
deficits and radii. Our approach avoids this degeneracy; however, 
it remains difficult to statistically determine whether scouring has 
taken place at large radii where the enclosed mass is much larger 
than Mbh ■ 

That being said, we do find that the median trend of mass 
deficit with radius, and absolute value inferred as a function of 
mass, can be approximately reproduced with the simple assumption 
that the "progenitor" profile represents a power-law continuation of 
the nuclear slope with some minimum dln/^/dlnr ~ —0.6, at radii 
where the slope flattens below this value. This is, of course, an as- 
sumption of the sort we are trying to avoid, but it does allow for 
estimation of the mass "deficit" on an individual galaxy basis, and 
our comparison here provides at least some justification for such a 
m ethodology. This als o exp lains why our results are si milar to those 
in lLauer et all ( l2007bl) and lMilosavlievic et alj ( I2OO2I) . It would be 
of considerable physical interest to see whether or not there is any 
physical motivation for such a "characteristic" slope - for example, 
it has been suggested that self-regulating inflow processes might 
set a characteristic slope around this value, in inflows dr iven by 
gravitational instabilities ( iHopkins & Ouataerlll2009l,l201(]|) . 

We stress that we agree with essentially all of the qualita- 
tive conclusions of the works discussed above. However, the re- 
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suiting quantitative distinctions and possible systematics are im- 
portant as tests of the theoretical models for core creation via BH 
binary mergers; very large AM* and/or apparent core radii ~kpc 
are actually quite difficult for models of core scouring to explain. 
Our more conservative study here therefore presents a natural and 
important complement to these works. We robustly set limits on 
the median scouring mass as a function of galaxy mass and radius 
without reference to any specific functional form for profile fitting. 
Moreover, to the extent that these measurements can be improved, 
they hold the potential to constrain not just the total mass that has 
been affected by scouring, but actually the scouring "kernel" - i.e. 
the magnitude of the effect as a function of radius, which can break 
a number of degeneracies between t heoretical models (see the dis- 
cussion in lPerets & Alexandeill200^ . 

Clearly, the very simple methodology presented here has sig- 
nificant limitations. We find even with a sample of ~ 200 well- 
observed spheroids that the statistical significance of our results is 
limited by the number of observed systems, because we can only 
consider systems in a narrow stellar mass range, and must in that 
range have a significant population of both cusp and core systems. 
As such, improving these constraints will require extending these 
samples with more observations of rare low-mass core or high-mass 
cusp ellipticals. Increasing the number of systems will also extend 
the radial range over which constraints can be applied; at some ra- 
dius R with an average enclosed mass M, (< R), the expected frac- 
tional offset between populations will be ~ Mbh/M, (< Ri) - this 
must be > 1/ \fN for a significant detection of the residual effects 
of "scouring" at these radii. In particular, determining whether or 
not the scoured mass is really ejected to large radii (or unbound) 
or simply kicked to moderate ~ 0. 1 — 1 kpc radii will require ex- 
tended samples along these lines. With such observations, study of 
the "scouring profile" (AM*(< R)) can be used to determine the 
differential contribution to scouring/BH binary coalescence from 
stars at different radii. This will allow direct observational probes 
of the loss cone, in particular of which stellar radii and orbits con- 
tribute. 

The methodology here also only constrains the median mass 
difference between populations - it cannot be used on an object-by- 
object basis. But with significantly larger samples, such that not just 
the median but the scatter and shape of each distribution M, (< R) 
can be determined, then the variation in scouring effects can be 
constrained (at least in a statistical sense). 

Independent tests are also important. For example, scour- 
ing is expected to preferentially eli minate stars on radial orbit s 
and leave a bia s for t angential orbits dOuinlan & HemquislllT997h . 
iGebhardt et J] ^20031) see tentative evidence for this in a limited 
sample of ellipticals; the major-axis radii within which the effect 
appears are generally ~ 0.5 — 3i?BH, similar to our estimates here. 
Ultimately, a combination of these independent constraints and the 
non-parametric constraints presented here can be used to deter- 
mined whether or not, and if so which, parametric approach(es) 
are most faithful to the true scoured mass, and motivate their use 
on an object-by-object basis. 
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